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We show that excitons in coupled quantum dots are ideal 
candidates for reliable preparation of entangled states in solid- 
state systems. An optically controlled exciton transfer process 
is shown to lead to the generation of Bell and GHZ states in 
systems comprising two and three coupled dots, respectively. 
The strength and duration of selective light-pulses for pro- 
ducing maximally entangled states are identified by both an- 
alytic, and full numerical, solution of the quantum dynamical 
equations. Experimental requirements to build such entan- 
gled states are discussed. 
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Quantum information, quantum computation, quan- 
tum cryptography and quantum teleportation represent 
exciting new arenas which exploit intrinsic quantum me- 
chanical correlations .EJ A fundamental requirement for 
the experimental realization of such proposals is the suc- 
cessful generation of highly entangled quantum states. 
In particular, coherent evolution of two quantum bits 
(qbits) in an entangled state of the Bell type is fun- 
damental to both quantum cryptography and quantum 
teleportation. Maximally entangled states of three qbits, 
such as the so-called Greenberger-Horne-Zeilinger (GHZ) 
statesa, are not only of intrinsic interest but are also of 
great practical importance in such proposals. New sys- 
tems and methods for the preparation and measurement 
of such maximally entangled states are therefore being 
sought intensively. Most of the theoretical and experi- 
mental activity to date hasJpeen associated with atomic 
and quantum-optic systemsLl'B. 

Solid-state realizations of such quantum-based phe- 
nomena have received little attention despite the fact 
that semiconductor nanostructures such as quantum dots 
(QDs), with quantum- mechanical electron confinement 
in all three directions, have been fabricated and studied 
by many groups.u Lrt^ddition, recent experimental work 
by Bonadeo et alE^B suggests that optically-generated 
electron-hole pairs (excitons) in semiconductor QDs rep- 
resent ideal candidates for achieving coherent wavefunc- 
tion control on the nanometer and femtosecond scales. 

In this paper we give a detailed prescription for produc- 
ing such entangled states in semiconductor quantum dot 
systems. We show that the resonant transfer interaction 
between spatially separated excitons can be exploited to 
produce such entanglement, starting from suitably ini- 
tialized states. The system requirements are realizable 
in current experiments employing ultrafast optical spec- 



troscopy of quantum dots. 

When two quantum dots are sufficiently close, there 
is a resonant energy-transfer process originating from 
the CoulomhL-interaction whereby an exciton can hop 
between dotsu. Experimental evidence of such energK- 
transfers between quantum dots was reported recentlyQ; 
the resonant process also plays a fundamental role in bi- 
ological and orsapic systems, and is commonly called the 
Forster processll^. Unlike usual single-particle transport 
measurements, the Forster process does not require the 
physical transfer of the electron and hole, just their en- 
ergy. Hence it is relatively insensitive to effects of im- 
purities which lie between the dots. Consider a system 
of A'' {N = 2, 3 or 4) identical, equispaced QDs contain- 
ing no net charge, radiated by long-wavelength classical 
light. Ignoring any constant energy terms, the following 
Hamiltonian describes the formation of single excitons 
within the individual QDs and their inter-dot hopping: 
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Here ejj (/ijj) is the electron (hole) creation operator in 
the n'th QD. The QD band gap is e while V represents 
the interdot Coulomb interaction and hence the Forster 
process. The dots are equidistant from each other, i.e. 
N — 2 dots on a line, = 3 dots at the vertices of 
an equilateral triangle, = 4 dots at the vertices of a 
pyramid; hence V is not a function of n or n'. The time- 
dependence of E{t) describes the pulse shape, while the 
magnitude includes the electron-photon coupling and the 
incident electric field strength. The Hamiltonian can be 
easily manipulated using quasi-spin operators 
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These satisfy the usual angular momentum commuta- 
tion relationships: [J+, J_] = 2Jz, [J±tJz] = T<^±- The 
Hamiltonian now takes the form 



H{t) = eJ, - V{J^ ~ Jl) + E{t)J+ + E*{t)J- 
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II{t) contains a non- linear term that can be exploited to 
generate entangled states. 
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We consider a rectangular radiation pulse, starting at 
time t = with central frequency to, given by E{t) = 
Acos{ujt): such a pulse is relatively straightforward to 
achieve experimentally. The time evolution of any initial 
state under the action of H in Eq.(3) is easily perfomjed. 
by means of the pseudo 1/2-spin operator formalismU'E£l 
Single transition operators are defined by 
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Here / denotes the identity matrix in the subspacc J — 1. 

In the absence of light, the energy levels of the system 
are given hy E = -Auj-V, -2V and Aw - V. Note that 
the energy separation between states |0) and |2) is unaf- 
fected by the inter-dot interaction V. Now we consider 
the action of a pulse of light at resonance, i.e. Acu = 0, 
and amplitude given hy A <^ V. We assume that the 
decoherence processes are negligibly small on the time 
scale of the evolution (see later). The density matrix at 
time t becomes 



where r — s denotes the transition between states |r) 
and \s) within a given J subspace. The three opera- 
tors belonging to one particular transition r — s obey 
standard angular momentum commutation relationships 
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where (a, (3, 7) represents a cyclic which exhibits the generation of coherence between vac- 



permutation of ix,y,z). Note that operators belonging ^um and bicxciton states through the operator J, 
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to non-connected transitions commute: 

with a, P — x,y or z. In order to gain physical insight 
into the dynamics of such a multi-exciton problem, some 
approximations are necessary: a common assumption, 
valid when e ^ y, is the so-called rotating wave ap- 
proximation U — e^*'^'^^*. Suppressing rapidly oscillat- 
ing terms, the Hamiltonian in the rotating frame (RF) 
becomes 

Hr = U^HU + i—-U = AujJ^ - AJ^ ~ V(J^ - J^) (5) 
at 

where Aw = e — uj denotes the off-resonance condition. 
We now show that this Hamiltonian leads to the genera- 
tion of entangled states from suitably initialized states. 

Two coupled QDs: Bell states. Here we describe the 
light excitation procedure to obtain a maximally entan- 
gled Bell state of the form \^) = |00) + e""^]!!) with 
(1) denoting a zero-exciton (single exciton) QD. The 
phase angle (j) can be arbitrary. In order to highlight 
the physical aspects of the procedure, we first derive an 
approximate analytical solution of the dynamical equa- 
tion governing the system's matrix density. Starting 
with the initial condition representing the vacuum of 
excitons, |J = = —1), only the J = 1 subspace 

is optically active while the J — subspace remains 
dark. We choose the basis of eigenstates of and 
J„ {|0) = |J = 1,M = -1), |1) = |J = 1,M = 0), 
|2) = |J = l,il/ = 1)}, as an appropriate representation 
for this problem. |0) represents the vacuum for exci- 
tons, |1) denotes a symmetric delocalized single-exciton 
state while |2) represents the biexciton state. The RF- 
Hamiltonian and initial density matrix can be expressed 
in terms of pseudo-spin operators as follows: 
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which oscillates at frequency 002 = A'^ /V. 

The Bell state l^-) = |00) + e*'^|ll) has a correspond- 
ing density matrix pBeii = 1/3 + — + 
cos((/)) J°~^ — sin(0) J°~^. Comparing this last equation 
with Eq.(8), we see that the system's quantum state at 
time T2 — TrV/2A'^ corresponds to a maximally entan- 
gled Bell state of the desired form with cj) = 7r/2. The 
time evolutions of populations and coherences for an ini- 
tial vacuum state are plotted in Fig. 1. The evolu- 
tion of populations of the vacuum poo and the biexciton 
P22 states are shown in Fig. la. Clearly our approxi- 
mate analytic calculation describes the system's evolu- 
tion very well when compared with the exact numer- 
ical solution (Fig. la). Figure lb shows the overlap, 
0{t) = Tr [pBeUp{t)], between the maximally entangled 
Bell state and the one obtained by applying a rectangu- 
lar pulse of light at resonance. The thick solid line (Fig. 
lb) describes 0{t) with a maximally entangled Bell state 
in the rotating frame, while the thin solid line (Fig. lb) 
represents the overlap with a Bell state transformed to 
the laboratory frame: obviously the RF case corresponds 
to the amplitude evolution of the laboratory frame sig- 
nal. The dashed line illustrates the approximate solution 
overlap in the RF. The approximate solution works very 
well, supporting the idea that a selective Bell pulse of 
length T2 = nV/2A'^ can be used to create a Bell state 
(0 = n/2) in the system of two coupled QDs. The same 
conclusion can also be drawn from the time evolution of 
the overlap between the exact Bell-state density matrix 
and the one obtained directly from the numerical calcu- 
lation. Therefore, the existence of a selective Bell pulse 
is numerically confirmed. 

Three coupled QDs: GHZ states. Next consider three 
quantum dots of equal size placed at the corners of an 
equilateral triangle. We can consider the J = 3/2 sub- 
space as being the only optically active subspace: the 
other two J = 1/2 subspaces remain optically dark. We 
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work in the basis set |J = 3/2, M), {|0) = |3/2,-3/2), 
|1) = 13/2,-1/2), |2) = 13/2,1/2), |3) = |3/2,3/2)}, 
where |0) is the vacuum state, |1) is the symmetric de- 
locahzed single-exciton state, |2) is the symmetric delo- 
cahzed biexciton state and |3) is the triexciton state. In 
the absence of radiation the energy levels are —3{Auj + 
V)/2,-{Au; + 7V)/2, (Aw - 7V)/2 and 3(Aw - V)/2. In 
terms of pseudo-spin operators, the Hamiltonian in RF, 
including the radiation term, is now given by 



-Aw(3Jr3 + J^-^) + 2V{J°,-^ - 2J2-3) 



(9) 



Two kinds q£ maximally entangled GHZ states have to 
be consideredtj: 

(i) the entangled state between vacuum and triexciton 
states given by \GHZ)i = -ij(|0) + e^'^|3)) or in terms of 

its associated density matrix by pci — I/^ + J°~^/2 — 
Jl~^/2 + cos((/)) J°~3 + sin(</)) J°-3, where / denotes the 
identity matrix in the J = 3/2 subspace. We now show 
that this state can be generated after an appropriate 7r/2- 
pulse. Starting with a zero-exciton state |0), the evolved 
state under the action of Hamiltonian Eq.(9) at reso- 
nance, i.e. Aw = 0, can be obtained in a straighforward 
way in the limit A/V <C 1 using the properties of pseudo- 
spin operators: 
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with uj^ — — dj^ + A and 

Clearly \GHZ)i (with 4> — 7r/2) can be generated with 
a 7r/2-pulse of length = ^-kV"^ /iA^. In Fig. 2a we 
show the overlap between the exact density matrix and 
that corresponding to state \GHZ)i. The dashed line 
shows the overlap using our approximate density matrix, 
Eq.(lO). 

(ii) the entangled state between a single exciton |1) and 
the biexciton |2) given by \GHZ)2 = ^(|1) -f e*'^|2)), or 
in terms of the corresponding density matrix pQ2 =1/4— 
J°- 72 + J2-3/2 -I- cos(</>) _^ sin(0) ji-2. In order to 
generate \GHZ)2 the initial condition must be modified. 
Using a suitably designed preparation sequence of pulses 
a new initial state, corresponding to a single exciton state 
|1), can be generated. Evolution of this new initial state 
under Hr (Eq.(9)) with Acu = 0, generates a density 
matrix at time t given by 



p(t) = i/-ij^i 
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where uj'^ = d 



d^ — A. Comparing this last result with 
the density matrix corresponding to a \GHZ)2 state we 



see that a 7r/2-pulse can be identified with a duration 
given by Tg = n/AA. Figure 2b shows the overlap be- 
tween p{t) (Eq.(ll)) and pG2- There is good agreement 
between our approximate result, Eq.(ll), and the exact 
one. 

We emphasize that the two maximally entangled GHZ 
states considered above have very different frequencies. 
In the limit A/V <C 1 the state \GHZ)i oscillates at 
the frequency W3 - 3A^/8V'^ while the state \GHZ)2 
oscillates at the larger frequency ~ 2A. This fea- 
ture should enable each of these maximally entangled 
GHZ states to be manipulated separately in actual ex- 
periments, even if the initial state is mixed. Furthermore 
we note that, after the preparation step, the system is 
evolving under the action of the Hamiltonian in Eq. (5) 
with Auj = A = 0: each one of the maximally entan- 
gled states discussed in this work are eigenstates of this 
remaining Hamiltonian. Hence in the laboratory frame 
\Bell) oscillates at frequency 2e, \GHZ)i oscillates at 
frequency 3e and \GHZ)2 oscillates at frequency e. 

Experimental observation of these Bell and GHZ states 
should be possible with present ultrafast semiconduc- 
tor optical techniquesQEl. For instance, for GaAs QDs 
e = 1.4eV which implies a resonance optical frequency 
u! = 2xl0"'^^sec~"'^. From the results above, it follows 
that in order to generate maximally entangled exciton 
states, 7r/2-pulses with sub- picosecond duration should 
be used. Femtosecond optical spectroscopy iflr-|GaAs- 
based nanostructures is now currently in usall'El. On 
the other hand, wide-gap semiconductor QDs, like ZnSe- 
based QDs, should do better because of the shorter re- 
quired 7r/2-pulse length: e = 2.8eV which leads to an 
optical resonance frequency of w = 4xl0^^sec^^. Fem- 
tosecond-spectroscopy is also currently available for this 
system.Ej In addition, the corresponding increase in the 
effective gap will yield a larger exciton binding energy: 
typical decoherence mechanisms (e.g. acoustic phonon 
scattering) will hence become less effective. A surprising 
conclusion of our results is that entangled-state prepa- 
ration is facilitated by weak light fields (i.e. A <^ V): 
strong fields cause excessive oscillatory behavior in the 
density matrix. This paper has considered the relatively 
straightforward experimental situation of global excita- 
tion pulses, i.e. pulses acting simultaneously on the en- 
tire QD system. However, otherj^ossibilities exist such 
as near-field optical spectroscopylij which allows the op- 
tical excitation and detection of individual QD signals. 
In this way, maximally entangled states with different 
symmetries can be obtained. 

Finally there is the well-known but difficult problem of 
decoherence. Existing solid-state proposals for quantum 
computers include quantum gates in coupled QDs based 
on electron spin effectsEJ and electronic charge effectstZl. 
Although electronic charge effects are subject to phonon 
decoherence, subdecoherent information may still be en- 
coded in such a quantum-dot array as described recently 
by Zanardi et alJia The present paper has focussed on 
exciton-based systems, since these have already been 
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shown to exhibit gOjOd coherence properties up to the pi- 
cosecond time-scaletl. Phonon decoherence wiU therefore 
be relatively unimportant on this time-scale. In addition, 
since no inter-dot transport of particles occurs, scatter- 
ing due to impurities etc. lying between the dots will be 
negligible. A detailed analysis of all possible decoherence 
times is beyond the scope of the present work, but will 
be addressed elsewhere. 

In summary, we have shown how maximally entangled 
Bell and GHZ states can be generated using the optically 
driven resonant transfer of excitons between quantum 
dots. Selective Bell and GHZ pulses have been identi- 
fied by an approximate, yet accurate, analytical approach 
which should prove a useful tool when designing exper- 
iments. Exact numerical calculations confirm the exis- 
tence of such 7r/2-pulses for the generation of maximally 
entangled states in coupled dot systems. 
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Figure Captions 

Figure 1: (a) Population of the vacuum state poo Siiid 
biexciton state P22 in two coupled QDs, as a function 
of time, (b) Time-evolution of overlap with maximally 
entangled BeU state, e = 1, V = e/10 and A = V/5. 
Thin solid line shows exact numerical result in the labo- 
ratory frame. Thick solid line in (b) represents the exact 
numerical solution in the rotating frame. Dashed line 
shows approximate analytic result. 

Figure 2: Time-evolution of overlap with maximally en- 
tangled GHZ states, (a) \GHZ)i and (b) \GHZ)2 under 
the action of a rectangular pulse of light at resonance. 
e = I, V = e/10 and A = 2V/5. Thick sohd fine rep- 
resents exact numerical solution. Dashed line shows ap- 
proximate analytic result. 
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